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ABSTRACT 

The metal deficient (Z = Z /41) Blue Compact Dwarf Galaxy (BCD) SBS 0335-052 
was observed with ISOCAM between 5 and 17/um. With a L^m/L^ ratio of 2.15, 
the galaxy is unexpectedly bright in the mid-infrared for such a low-metallicity object. 
The mid-infrared spectrum shows no sign of the Unidentified Infrared Bands, which we 
interpret as an effect of the destruction of their carriers by the very high UV energy 
density in SBS 0335-052. The spectral energy distribution (SED) is dominated by a 
very strong continuum which makes the ionic lines of [Siv] and [Nein] very weak. 
From 5 to 17 /um the SED can be fitted with a grey-body spectrum, modified by an 
extinction law similar to that observed toward the Galactic Center, with an optical 
depth of Ay ~ 19-21 mag. Such a large optical depth implies that a large fraction 
(as much as ~ 75%) of the current star-formation activity in SBS 0335-052 is hidden 
by dust with a mass between 3xlO 3 M and 5xlO 5 M . Silicate grains are present 
as silicate extinction bands at 9.7 and 18 /xm can account for the unusual shape of 
the MIR spectrum of SBS 0335-052. It is remarkable that such a nearly primordial 
environment contains as much dust as galaxies which are 10 times more metal-rich. 
If the hidden star formation in SBS 0335-052 is typical of young galaxies at high 
redshifts, then the cosmic star formation rate derived from UV/optical fluxes would be 
underestimated . 



Subject headings: stars: formation - dust, extinction - galaxies: compact - galaxies: 
individual (SBS 0335-052) - infrared: ISM: continuum - infrared: ISM: lines and bands 



1 Based on data obtained with ISO, an ESA project with instruments funded by the ESA member states (especially 
the PI countries: France, Germany, the Netherlands, and the United Kingdom) with the participation of ISAS and 
NASA. 
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1. Introduction 

Galaxy formation is one of the most fundamental unsolved problems in astrophysics. Much 
effort has gone into the search for primeval galaxies (hereafter PGs), i.e. galaxies undergoing their 
first major burst of star formation, at redshifts larger than ~2. Several objects have been put 
forward as possible PGs, ranging from high-redshift radio-galaxies to Lya emitters (e.g. Steidel 
et al. 1996; Yee et al. 1996). However, most of these candidate PGs appear to already contain 
a substantial amount of heavy elements (as indicated for example by the presence of P Cygni 
profiles), implying previous star formation and metal enrichment. 

SBS 0335-052 (a(1950) = 03^ 35 m 15.2 s , 5(1950) = -05° 12' 25.9") is a relatively nearby Blue 
Compact Dwarf Galaxy (hereafter BCD) with an absolute magnitude M#= -16.7 and which 
appears to be undergoing its first burst of star formation. With a metallicity of only Zq/41 
(Izotov et al. 1997), it is the second most metal-deficient galaxy known, after I Zw 18 (Z0/5O). 
With a redshift of 0.0136 and a Hubble constant of 75 km s _1 Mpc -1 , it is at a distance of 54.3 
Mpc (1" = 263 pc). Thuan, Izotov, & Lipovetsky (1997) and Papaderos et al. (1998) have found 



that the stars in SBS 0335-052 are not older than ~ 100 Myr, making it a truly young galaxy. 



Thuan fc Izotov (1997)| suggest that the large HI envelope within which SBS 0335-052 is embedded 
may be truly primordial. Yet, despite its youth and extremely low metallicity, HST images of 
SBS 0335-052 clearly show dust patches mixed in with the six super-star clusters where most of 
the star formation is occurring (Thuan et al. 1997). The presence of dust, in combination with 
the intense ultraviolet radiation field from the many young stars in SBS 0335-052 suggests that 
there may be detectable Mid-Infrared (MIR) emission, as the dust will reprocess the UV starlight 
and reemit it in the infrared. We have therefore obtained MIR observations of SBS 0335-052 with 
the Infrared Space Observatory (ISO) to study the properties of dust and star formation in a truly 
metal-deficient environment, similar to those prevailing at the epoch of galaxy formation. 



2. Observations and Data Reduction 

The observations were obtained with ISOCAM (Cesarsky et al., 1996a), the MIR imaging 
instrument aboard ISO (Kessler et al., 1996). They consist of two different sets: a set of 
broad-band filter maps, and a set obtained by imaging spectroscopy. 

The broad-band maps were obtained with a spatial sampling of 3" per pixel and in a 3x3 
raster mode with a 3-pixel displacement in each direction, giving a total field of view of 114" x 114". 
The raster mode increases the sensitivity and improves the flat-field correction since in the center 
of the map, each sky position is imaged by 9 different camera pixels. The maps were obtained 
with the LW9, LW10, LW8, LW6 and LW2 filters, centered respectively at 14.9, 12.0, 11.3, 7.7 and 
6.7 fiia (Table |l|). Note that the LW10 filter is identical to the IRAS 12 /mi filter. The integration 
time was 10.08 s per readout for all filters, except LW10 for which an integration time of 5.04 s was 
used. Total on-source times were 860 s for LW10 and 1540 s for the other filters. 



- 3 - 



Imaging spectroscopy was performed using the Circular Variable Filter (CVF) facility of 
ISOCAM. This mode produces images of the full ISOCAM field of view (192" x 192") in the 
2-17 yum wavelength region, with a sampling of 6" per pixel and a spectral resolution A/AA of ~ 40. 
Due to the faintness of the source we obtained observations only in the 8.78-17.34 fim wavelength 
range. The individual integration time was 5.04 s and the total on-source time was 4800 s. 



The various steps of ISOCAM data reduction are detailed in Stark et al. (1998) and we will 



only describe here aspects that depart from the methods discussed in that paper: dark current 
subtraction and transient correction. 

The dark current is known to show secular trends which depend on the position of the satellite 
in its orbit and on the time elapsed since launch. These trends can be accurately modelled (see 
Biviano et al. 1998), allowing nominal dark subtraction. 

Transient corrections using the inversion algorithm of Abergel et al. (1996) were 
applied successfully to the broad-band observations. For the CVF observations, the known 
oversimplification of the method combined with the intrinsic faintness of the source results in 
errors of similar magnitude for both corrected and uncorrected CVF data. The main effect of 
the correction is to add an artifact at the start of the spectrum (around 17 /im) while leaving 
the rest of the scan unchanged. In particular, synthesis of the LW8 and LW9 flux densities from 
corrected or uncorrected spectra gives flux densities 1.3 times larger than the observed broad-band 
flux densities. That the CVF produces higher flux densities than the broad-band filters is to be 
expected: the operational setup of ISOCAM is such that the scan starts with a short exposure 
through the LW2 filter with a scale of 6" per pixel, giving a much higher illumination on the 
detector than the scan itself. As the source is very faint, the detector would need a much longer 
time than our adopted exposure time to stabilize down to its true level. This is not the case for 
the broad-band maps as the broader filters allow many more photons to reach the detector, thus 
speeding the stabilization process. As a result, we have adopted for the remaining discussion the 
CVF spectrum uncorrected for transient effects and scaled down by a factor of 1.3 . 

Figure 1 displays the spectral energy distribution of SBS 0335-052. The flat-fielding and 
photometry of the CVF scan were made following Aussel et al. (1998). Conversion from camera 



units to mJy was performed using the calibration factors given in the ISOCAM Cookboo% HST 
images (Thuan et al. 1997) reveal that most of the star formation in SBS 0335-052 occurs within a 
region of ~2.5" in size, so that the individual super-star clusters in SBS 0335-052 are not resolved 
by ISOCAM. Thus the flux densities given in Table |l| and plotted in Figure 1 are integrated over 
the whole star-forming region. 



3. Discussion 



2 available at 



http : //isowww. estec. esa.nl/manuals/iso_cam/ 
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3.1. Global properties 

Since the ISOCAM LW10 filter is equivalent to the IRAS 12 ^urn filter, it is interesting to see 
how the integrated properties of SBS 0335-052 compare with those of other BCDs and irregular 
galaxies. We have thus compiled IRAS and Ha fluxes for the Thuan &: Martin (1981, hereafter 
TM81) catalog of BCDs, as well as for the sample of BCDs and irregular galaxies in Sauvage, 
Thuan, & Vigroux (1990) to perform the comparison. 

With a 12 fim to blue luminosity ratio L12/LB of 2.15=L0.06, SBS 0335-052 is comparable 
to the most MIR-bright objects in the TM81 catalog. In that sample the mean L^/Lb ratio 
(<logLi2/Le> =-0.12±0.27), while already twice as high as that for spiral galaxies in the CfA 
catalog (Sauvage & Thuan 1994), is a factor of 3 smaller than the ratio for SBS 0335-052. This 
shows that even though the BCD is ~ 10 times more metal-poor than the galaxies in TM81, it 
nevertheless contains a significant amount of dust heated by an intense UV radiation field. 

If we compare the 12 11m. luminosity to a tracer of star formation such as the Ha luminosity 
(Table |2[), SBS 0335-052 appears to be normal. Although total Ha fluxes are scarce for BCDs, the 
available data give <logLi2/L# a > = 0.96±0.70 for the above samples, as compared to 0.65±0.03 
for SBS 0335-052. Thus, even though the metallicity of the BCD is exceptionally low, and the 
galaxy unusually MIR-bright, its star-forming properties appear to be normal when compared to 
other BCDs. 



3.2. The mid-infrared spectrum 

When the MIR spectral energy distribution of SBS 0335-052 (Figure 1) is compared to MIR 
spectra of other star-forming galaxies (see e.g. Metcalfe et al. 1996; Vigroux 1997), two facts stand 
out : (1) there is no emission from the so-called Unidentified Infrared Bands (hereafter UIBs, 
usually attributed to carbon-based dust, see e.g. Papoular et al. 1996; Puget & Leger 1989). In 
particular the very strong band at 11.3 //m is completely absent from the spectrum; and (2) there 
are no evident fine structure ionic lines, even though neon as well as sulfur lines are usually quite 
bright in starburst galaxies. 

In Figure 1, a hint of the [Siv]A10.5/iin and [Nein]A15.6/xm lines can be seen. There is a feature 
near the position of the [Neii]A12.8/im line, but it is centered at 12.8 //m, while the redshifted line 
falls at 12.98 /im. The quality of our spectrum prevents us from actually measuring line fluxes. 
Instead we compute upper limits as twice the local rms noise times the instrumental profile width 
at 20% of peak intensity at the line location. We obtain upper limits of 5.6x10" ergs -1 cm -2 
for the [Siv] line, 4.5xl0 -14 ergs -1 cm -2 for the [Nell] line, and 5.4 xl0 -14 ergs -1 cm -2 for the 
[Nem] line. 



Using the photoionization models of Stasinska (1990) with a metallicity equal to 0.02 that 
of the sun, and with an integrated H/5 flux in SBS 0335-052 of 12.1 xlO -1 ergs -1 cm -2 (the flux 
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given by Izotov et al. 1997 was multiplied by 4 to account for aperture and extinction effects as 
suggested by Thuan & Izotov 1997), we predict the line intensities of [Siv], [Nein], and [Nell] to 
be 5.0 xlO -14 erg s -1 cm -2 , 1.6 xl0~ 14 ergs -1 cm -2 , and 7.7xl0~ 17 ergs -1 cm~ 2 respectively. The 
predicted intensities are thus completely consistent with the upper limits. We conclude, therefore, 
that the weakness (or absence) of ionic fine-structure lines in SBS 0335-052 when compared to 
other star-forming galaxies can be explained by a very strong continuum which decreases the 
equivalent widths of the lines, making them difficult to detect. 

None of the usual UIB features are seen, suggesting that their carbonaceous carriers are 
absent or have been destroyed. As the star-forming region, 660 pc in diameter, contains 5000 



equivalent 07 stars (Izotov et al. 1997), the average energy density is ^lOeVcm . [Puget & 



Leger (1989) suggest that at this radiation density, the abundance of the band-emitting dust is 
reduced by a factor of ~10. Therefore the absence of the bands can be understood as a destruction 
effect. This is probably enhanced by the very low metallicity of the galaxy which allows harder 
UV photons to travel further into the interstellar medium. 

The origin of the continuum emission is more problematic. We observe what could be 
interpreted as either a broad emission feature at around 14 /im, or a broad absorption feature at 
A ^ 16 /im. We favor the second hypothesis as we know of no candidate which would show such an 
emission feature at ~ 14 /im. Very small carbonaceous grains, usually thought to be responsible 
for the MIR continuum (Desert, Boulanger, & Puget 1990), produce featureless emission. Silicates, 
the other well-known component of dust, can produce emission features of various shapes, but 
these are centered around 10 /im and even in circumstellar disks never extend beyond 12 /im (see 
e.g. Sitko et al. 1998). On the other hand, silicate extinction bands at 9.7 and 18 //m can provide 
an explanation for the strange shape of the MIR spectrum of SBS 0335-052. 

Testing this hypothesis is not straightforward as the expected shape of the dust spectrum in 
this wavelength range is difficult to constrain. Apart from the stochastic heating of macromolecules 
which produce the UIBs, MIR emission is thought to be produced by very small grains in a 
thermal regime intermediate between equilibrium and stochastic heating (see e.g. Desert et al. 
1990). We have therefore attempted to fit the spectral energy distribution (SED) of SBS 0335-052 
with different screen models where the emission spectrum was successively: (1) a black-body of 
variable temperature modified by an emissivity law of the form i u aB u (T) i^ 1-2 ; (2) a power-law 
spectrum; and (3) the continuum MIR spectrum of the Galactic HII region M17 with UIBs and 
ionic lines removed (Cesarsky et al. 1996b). The free parameters of the underlying spectra include 
in all cases a scaling factor, in addition to the temperature in case (1) and the power-law slope in 
case (2). A further uncertainty comes from the extinction curve to be used. We have tried two 



different laws, the standard MIR extinction law derived by Draine (1989) , and that observed in 



the direction of the Galactic Center by Lutz et al. (1996")| The second law differs from the first by 



its much higher A\/Ay in the range 2-8 /im. This difference is probably due to the neglect of ice 
coating on grains in the praine (1989) model. The free parameter associated with the extinction 



curve is simply the absolute extinction at a reference wavelength. Our fitting procedure searches 
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for best fits to the CVF scan data only. There is no attempt to fit the broad-band data outside of 
the CVF range. However, a successful model will be one that also accurately predicts them. 

Figure 2 shows examples of the three types of fits. Both a power law and the M17 spectrum 
fail to reproduce the observed SED. Power laws modified by extinction overproduce emission 
shortward of 8 /xm (Figure 2a) while the M17 spectrum overestimates emission longward of 14 /im 
(Figure 2b). On the contrary, a black-body spectrum modified by an emissivity law proportional 
to v 1 ' 5 and extinguished by a screen of dust, gives an excellent fit to the SED of SBS 0335-052 
(Figure 2c). The goodness of the fit depends very little on the exponent of the emissivity law, so 
its value cannot be used to constrain the nature of the emitting dust. Instead, the fit is sensitive to 
the shape of the extinction law. It is impossible to obtain a good fit for the LW2 and LW6 bands 



with the standard extinction curve of Draine (1989) . The predicted spectrum always overproduce 



the emission in these bands by factors of ~2. On the contrary, excellent fits of these bands can 



be obtained with the Lutz et al. (1996) extinction curve. The temperature obtained for the 
black-body curve is in the range 240-260 K, but probably carries little physical meaning since the 
MIR-emitting grains are not likely to be in thermal equilibrium. 

The information we derive concerning the dust extinction is more physical. Given the 
uncertainties in our data, we obtain an equivalent Ay of 19-21 mag. Note that since the screen 
model maximizes Ay for a given column density of dust, the column density derived from our fit 
is likely a lower limit to the true column density if emitting and absorbing dust grains are mixed. 

Such a high value of extinction is in sharp contrast to the values (Ay < 0.6 mag) derived by 



Izotov et al. (1997) from optical spectrophotometric observations. The large difference between 
the optical and infrared extinctions implies that the heating sources responsible for the infrared 
emission detected by ISO are probably too deeply embedded in dust to be detected in the visible 
and that, as in the case of the Antennae galaxies NGC 4038/39 (Mirabel et al. 1998), most of 
the infrared emission is powered by invisible star clusters. Indeed, correcting the Ha flux for 
20 mag of visible extinction would lead to unrealistically high star formation rates which would 
be incompatible with the non-detection of SBS 0335-052 at 1.4 GHz by the NVSS (Condon et al. 
1998), or in the far-infrared by IRAS. 

The LW10 luminosity of 2.5 xlO 8 L is equivalent to the bolometric luminosity of 2500 07 
stars, i.e. half the equivalent number of stars in the visible super-star clusters. Given that in 
the BCD sample of TM81, <logLi 2 /L F i R > = -0.85±0.09, it is probable that in fact some 18000 
equivalent 07 stars are required to power the total infrared emission. This implies that even in 
very low-met allicity environments such as those characterizing PGs, a significant fraction, i.e. ~ 
3/4, of the total star formation activity of a galaxy can be effectively hidden from UV/optical 
observations. 

From the extinction derived in the MIR, we can estimate the dust mass. As MIR extinction 
is essentially due to silicates (Draine &: Lee 1984), this result will likely be a lower limit as some 
carbon-based dust (such as the grains responsible for the continuum) can be present and still 
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contribute little extinction in the MIR. Using the dust model of Draine fc Lee (1984) , we obtain 
a dust surface density of 1.5 M pc~ 2 . To derive a dust mass, we need to estimate lower and 
upper bounds to the dust spatial extent. For a lower bound, we assume that the dust is only 
associated with hidden super-star clusters (SSCs). O'Connell, Gallagher, & Hunter (1994) found 
that SSCs typically have a very compact core with a half-light radius of ~3pc, embedded in 
considerably more extended halos with diameters ~30pc. If we adopt the latter value, and 
consider that at least 3 SSCs are hidden (they have the power of 2500 equivalent 07 stars needed 
to power the LW10 luminosity), we derive a dust mass of ~3.2 x 10 3 M . This is very probably 
an underestimate of the true dust mass as there are likely many more hidden SSCs. Furthermore 
dust is probably not just associated with the SSCs, but rather mixed throughout the star-forming 
region, especially if it has a supernova origin as argued later on. 

A likely upper limit to the spatial extent of the dust is the size of the region where reddening is 
observed by HST (Thuan et al. 1997), i.e. 660 pc. This gives a dust mass of 5xl0 5 M . With the 
Hi mass given in Table ||, the gas-to-dust mass ratio of SBS 0335-052 is then between ~ 2 x 10 3 
and ~ 3 x 10 5 , i.e. much higher than in the Galaxy, which is not too surprising for such a low 
metallicity object (Lisenfeld & Ferrara 1998). More extraordinary is the fact that the dust mass 
in SBS 0335-052 can be of the same order of magnitude as that in BCDs which are on average ~ 
10 times more metal-rich. For comparison <logM^ US f> = 4.4±0.6 M in the TM81 BCD sample. 

As discussed by [Thuan et al. (19971 , because stars in SBS 0335-052 are not older than ~100 
Myr, there is not enough time for the silicate dust to be made in the envelopes of red giant stars. 
Rather silicate grains probably condensed out of the numerous supernova ejecta present in the 
BCD (see e.g. [Lucy et al. (1991)| on SN1987A; Dwek et al. 1992; Wooden et al. 1993). We can 



check for the plausibility of this hypothesis by adopting, for example, the supernova silicate dust 
mass input rate of 0.5 M pc -2 Gyr -1 obtained by pwek (1998)| for the Galactic Center at 0.1 Gyr 
,the age of SBS 0335-052. Assuming the diameter of the dust-forming region to be 660 pc, we 
derive a silicate dust mass of ~2xl0 4 M , in the range of the above estimates. 



4. Conclusions and Implications 

Although SBS 0335-052 is one of the most metal-deficient galaxies known (Z /41), it is 
unexpectedly bright in the MIR range, implying a silicate dust mass in the range 10 3 -10 5 M . The 
MIR spectrum shows no sign of UIB carriers which are probably destroyed. Despite a difference 
of a factor of 40 in metallicity, the Galactic Center MIR extinction produces the best fit to the 
spectrum. This is quite different from the situation in the UV where the extinction law shows 
a strong dependence on metallicity (Fitzpatrick 1989). This difference is probably due to the 
fact that the abundance and spectral properties of the carbon-based dust responsible for the UV 
extinction and UIBs are more metallicity-dependent than those of the silicate grains responsible 
for the extinction in the MIR. A possible explanation is that the silicate grains are more resistant 
to photo-destruction than carbon-based dust. 
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The derived extinction is quite high, Ay ~ 20 mag. Given that the total MIR luminosity 
already requires the bolometric luminosity of 50% more young stars than are seen in the galaxy, the 
total star formation rate (SFR) as derived from the optical or UV luminosities must underestimate 
the true star formation rate by at least 50%, and more likely by a factor of 4 as argued above. 
Thus if SBS 0335-052 is a representative example of PGs, the cosmic star formation rate will be 
systematically underestimated if based only on UV and optical fluxes. This is in fact the result 
obtained by Flores et al. (1998) in their ISO survey of distant galaxies. They found that the 
cosmic SFR derived from FIR luminosities is ~ 2 to 3 times higher than the SFR estimated 
previously from UV/optical fluxes (Madau et al. 1996). 



The ISOCAM data was analyzed using the software package "CIA", a joint development by 
the ESA Astrophysics Division and the ISOCAM consortium. The ISOCAM consortium is led by 
the ISOCAM PI, C. Cesarsky, Direction des Sciences de la Matiere, C.E.A, France. T.X.T. has 
been partially supported by NASA grant JPL961535. We acknowledge useful conversations with 
Yuri Izotov. 
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Fig. 1. — Spectral energy distribution of SBS 0335-052 in the mid-infrared. The crosses represent 
broad-band observations at 6.7, 7.7, 11.3, 12.0 and 14.9 /j,m. The horizontal bars show the width of 
the broad-band filters. The la errors bars are smaller than the crosses. Broad-band flux densities 
are derived assuming an intrinsic spectral shape f u oc v~ x and integrating over the point spread 
function. The thick line is the CVF spectrum without transient correction and scaled down by 
a factor 1.3 (see text for details). The dashed lines show ± la errors. The spectrum has been 
corrected for effects of the varying point spread function with wavelength. Expected positions of 
the [Siv]A10.5/im, [Nell]A12.8/xm, and [Nelll]A15.6/um emission lines are indicated. 
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Fig. 2. — Best-fit models (dashed lines) for the spectral energy distribution of SBS 0335-052. Three 
different types of emission spectra are considered: (a) a power-law of the form u ' 7 with an extinction 
Ay = 27 mag; (b) the continuum observed in M 17 (Cesarsky et al. 1996b) with Ay = 0; (c) a black- 
body of temperature 260 K, with an exponent of 1.5 for the emissivity law, and Ay = 21 mag. For 
all three fits the extinction curve of Lutz et al. (1996) gives the best results. To compare our fit 
to the broad-band data (crosses), we have integrated the fitted spectrum (dashed line) over the 
band-pass of the ISOCAM filters. The resulting flux densities are shown by open diamonds. Only 
model (c) is able to reproduce satisfactorily all broad-band observations. 
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Table 1. MIR flux densities of SBS 0335-052 



filter 


A a 


Ai/ a 


AA a 


Flux density 


la 


Luminosity 13 




fim 


10 12 Hz 


/im 


mJy 


mJy 


L 


LW9 


14.9 


2.29 


2.0 


35.2 


0.9 


7.4 10 7 


LW10 


12 


13.20 


7.0 


20.5 


0.6 


2.5 10 8 


LW8 


11.3 


2.71 


1.3 


20.7 


0.8 


5.2 10 7 


LW6 


7.7 


7.06 


1.5 


8.9 


0.5 


5.8 10 7 


LW2 


6.7 


16.18 


3.5 


5.5 


0.4 


8.2 10 7 



a Assuming a spectral shape f u <xv , as in the ISOCAM cookbook. 
b Assuming a distance of 54.3 Mpc. 



Table 2. Relevant data for SBS 0335-052 



Parameter Value 

L Ha a 5.61O 7 L 

Lff,j a 2.21O 7 L 

L B b 1.21O 8 L 

L B C 7.210 8 L® 

M ffi d 9.51O 8 M 

a Flux in the 1" slit of Izotov et al. 1997, corrected for extinction using Ay = 0.6 mag, and 
multiplied by 2 to account for aperture effects as suggested by Thuan &; Izotov 1997. 

b Computed from the absolute blue magnitude in Thuan et al. 1997. 

c Here Lb is expressed in units of the solar blue luminosity, where Lq/L 1 ^ ~6.25 . 

d Thuan et al. 1998. 



